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Abstract 

In this paper we study thermoconvective instabilities appearing in a fluid within a 
cylindrical annulus heated laterally. As soon as a horizontal temperature gradient 
is applied a convective state appears. As the temperature gradient reaches a crit- 
ical value a stationary or oscillatory bifurcation may take place. The problem is 



modelled with a novel method which extends the one described in |T0[ . The Navier 
Stokes equations are solved in the primitive variable formulation, with appropriate 
boundary conditions for pressure. This is a low order formulation which in cylin- 
drical coordinates introduces lower order singularities. The problem is discretized 
with a Chebyshev collocation method easily implemented and its convergence has 
been checked. The results obtained are not only in very good agreement with those 
obtained in experiments, but also provide a deeper insight into important physical 
parameters developing the instability, which has not been reported before. 



1 



1 Introduction 



The problem of thermoconvective instabilities in fluid layers driven by a temperature 
gradient has become a classical subject in fluid mechanics 0- It is well known that 
two different effects are responsible for the onset of motion when the temperature 
difference becomes larger than a certain threshold: gravity and capillary forces. 
When both effects are taken into account the problem is called Benard-Marangoni 
(BM) convection [|]. 

To solve numerically hydrodynamical problems in the primitive variable for- 
mulation raises two questions, one is the handling of first order derivatives for pres- 
sure and the other is finding its boundary conditions ^ which are not intuitive. 
In thermoconvective problems pressure usually is avoided, for instance in Refs. ^ 
the method of potentials of velocity is used to eliminate it from the equations. This 
technique raises the order of the differential equations and additional boundary con- 
ditions may be required. This is particularly troublesome in cylindrical coordinates 
where high order derivatives cause awkward difficulties. In Ref. the primitive 
variables formulation is used although the use of a particular spectral method allows 



the removal of pressure. In Ref. |]T0[ the linear stability analysis of two convection 



problems are solved in the primitive variables formulation and taking appropriate 
boundary conditions for pressure. In this article very good accuracy is obtained 
both for cartesian and cylindrical coordinates. This method has been succesfully 
applied to describe experimental results in cylindrical containers In this paper 
we extend this method to study a BM problem in a cylinder heated laterally. 

The physical set up (Fig. 1) consists of a fluid filling up a container bounded 
by two concentric cylinders. The upper surface is open to the air and the fluid 
is heated laterally through the lateral walls which are conducting. This set up 
corresponds to the experiment reported in Ref. |]l2l, and also it is comparable to 



similar experiments in rectangular containers |jT3|, In [|T^] a similar problem is 
treated theoretically, however in this work rotation is considered and the gravity field 
is along the radial coordinate. In the problem we treat, as soon as a slight difference 
of temperature is imposed between the lateral walls, a stationary solution appears 
which is called basic state. When the temperature of the walls is modified the basic 
state can become unstable and bifurcate to different patterns, both stationary and 
oscillatory. When the bifurcation is stationary, the emerging pattern consists of a 
3D structure of rolls whose axes are parallel to the radial coordinate (radial rolls) 
and in the oscillatory case, the axes are tilted with respect the radial axis. These 
results coincide with those obtained experimentally |T2| , [13], . 

The basic state is calculated solving a set of nonlinear equations whereby 
the first contribution is obtained by a linear approach. To improve the solution we 
expand the corrections of the unknown fields in Chebyshev polynomials, and we 
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pose the equations at the Gauss-Lobatto collocation points [|I6|. The bifurcation 
thresholds are obtained through a generalized eigenvalue problem with the same 
Chebyshev collocation method. The convergence of the method is studied by com- 
paring different expansions. 

The organization of the article is as follows. In Sec. 2 the formulation of the 
problem with the equations and boundary conditions is explained. In the third sec- 
tion the basic state is calculated and the results obtained are discussed for different 
physical conditions. In the fourth section the linear stability of the basic solution is 
performed and instabilities are studied for different parameters. In the fifth section 
conclusions are detailed. 



2 Formulation of the problem 

The physical set up considered is shown in Fig. 1. A horizontal fluid layer of depth 
d {z coordinate) is in a container bounded by two concentric cylinders of radii a and 
a+6 (r coordinate). The bottom plate is rigid and the top is open to the atmosphere. 
The inner cylinder has a temperature Tmax whereas the outer one is at T^in and the 
environment is at Tq. We define AT = T^ax — and ATh = T^ax — ^min, which 
are the main two parameters controlling instabilities in this problem. The system 
evolves according to the momentum and mass balance equations and to the energy 
conservation principle. In the equations governing the system Ur, and Uz are the 
components of the velocity field u of the fluid, T the temperature, p the pressure, 
r the radio vector and t the time are denoted. The magnitudes are expressed in 
dimensionless form after rescaling in the following form: r' = r/rf, t' = K,t/(P, 
u' = duj K^j)' = (P"pj {^pqkv) , B = (T — To) / AT. Here k is the thermal diffusivity, v 
the kinematic viscosity of the liquid and po is the mean density at the environment 
temperature Tq. 

The governing dimensionless equations (the primes in the corresponding fields 
have been dropped) are the continuity equation, 

V-M = 0. (1) 

The energy balance equation, 

+ n • V0 = V^e, (2) 

The Navier-Stokes equations, 

+ V)M = Pr (-Vp + V\ + — ^^ej , (3) 
V apoAT ) 
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where the operators and fields are expressed in cyhndrical coordinates |Q and 
is the unit vector in the z direction. Here the Oberbeck-Bousinesq approximation 
has been used. It consist of considering only in the buoyant term, the following 
density dependence on temperature p = pq\V — a{T — Tq)], where a is the thermal 
expansion coefficient. The following dimensionless numbers have been introduced: 

= (4) 

K KV 

where g is the gravity constant. Pr is the Prandtl number which is assumed to have 
a large value and R the Rayleigh number, representative of the buoyancy effect. 

2.1 Boundary conditions 

We discuss now the boundary conditions (be). The top surface is fiat, which implies 
the following condition on the velocity, 

Uz = 0, on z = 1. (5) 

The variation of the surface tension with temperature is considered: cr (T) = ctq ~ 
7 (T — Tq) , where is the surface tension at temperature Tq, 7 is the constant 
rate of change of surface tension with temperature (7 is positive for most current 
liquids). This effect supplies the Marangoni conditions for the velocity fields which 
in dimensionless form are, 

M 

dzUr + MdrQ = 0, dzU^ H (9,^6 = 0, on z = 1. (6) 

Here M = 'yATd/ {nupo) is the Marangoni number. In our particular problem 
the Rayleigh and the Marangoni numbers are related in the same way as in the 
experiments in Ref. M = 9.2 ■ 10~^R/(f, so buoyancy effects are dominant. 
The remaining boundary conditions correspond to rigid walls and are expressed as 
follows, 

Ur = = Uz = 0, on z = 0, (7) 
Ur = = Uz = 0, on r = a* , r = a* + 6*, (8) 

where a* = a/d and 6* = 6/d. 

For temperature we consider the dimensionless form of Newton's law for heat 
exchange at the surface, 

dzQ = -Be, on z = l, (9) 
where B is the Biot number. At the bottom a linear profile is imposed, 

^ / r a\ ATh , , 



4 



while in the lateral walls conducting boundary conditions are considered 

e = 1, on r = a*, (11) 
e = (-1 + ^)^ + 1, on r = a* + 6*. (12) 

/^From this boundary conditions it is clear how not only AT defining the Rayleigh 
number is involved in the problem, but also AT^. Due to the fact that pressure is 
kept in the equations, additional boundary conditions are needed. They are obtained 
by the continuity equation at 2; = 1 and the normal component of the momentum 



equations on r = a* , r = a* + S* and z = 0, . 

V • M = 0, on z = l, (13) 

p^_l / dUr _^ ^ dUr_ ^ U^dUr _^ ^ dUr_ _ _ 

\ dt " dr r dcj) ^ dz r j 

, A Ur 2 du^ , 

= -TT + J - ^^T-, on r = a,r = a +5, (14) 

or d<p 

Pr-^ ( ^ + u — + + u — ^ - 

\ dt ^ dr r d(f) ^ dz J 

= + Au^-b + RT, on z = 0, (15) 
dz 

where b = (fg/{Kh'), and A = r'^d / dr{rd / dr) + r~'^d'^ / d(j)^ + d'^/dz^. 



3 Basic state 

As soon as the lateral walls are at different temperatures and a horizontal tem- 
perature gradient is set at the bottom, a stationary convective motion appears 
in the fluid. In contrast to the classical Benard-Marangoni problem with uni- 
form heating, here the basic state is not conductive, but convective. In order to 
calculate it, for computational convenience, the following change has been per- 
formed: r' = 2r/6* — 2a*/ 6* — 1 and z' = 2z — 1, which transforms the domain 
rii = [a*, a* + 6*] X [0,27r] x [0,1] into n = [-1,1] x [0, 27r] x [-1,1]. After these 
changes, and since the basic state has radial symmetry (there is no dependence on 
(j)) the steady state equations in the inflnite Prandtl number approach become (the 
primes have been dropped), 

AdrP+C^Ur = A*Ur, (16) 

2d,p + h-RQ = A*u„ (17) 

GUr + AdrUr + 2d^u^ = 0, (18) 

UrAdre + 2uAe = A*e, (19) 
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where A* = A'^d^.+GAdr + Adl b = d^g/nu, A = 2d/5 and G (r) = 2d/{2a + 6 + r5). 
The boundary conditions are now as follows, 

u, = 2d,Ur + MAdrT = 23,0 + 39 = 0, on z = 1, (20) 

Gur + AdrUr + 2dzUz = 0, on z = 1, (21) 

= = 0, = (1 - r)/2, on z= -1, (22) 

2dzP + h- RQ = /\*u,, onz = -l, (23) 

= = 0, = 0, on r = 1, (24) 

AdrP + G'^Ur = A*Ur, on r = 1, (25) 

Uj. = Uz = 0, Q = 1, on r = —1, (26) 

AdrP + CV = ^*Ur, on r = -1. (27) 

3.1 Numerical method 

We have solved numerically Eqs. (|TH)-(|Ty|) together with the boundary conditions 



([20|)-(p7D by a Chebyshev collocation method. This approximation is given by four 
perturbation fields Ur (r, z) , Uz (r, z) , p (r, z) and (r, z) which are expanded in a 
truncated series of orthonormal Chebyshev polynomials: 

N M 

Ur (r, z) = (^nmTn (r) {z) , (28) 

n=0 m=0 
N M 

Uz (r, z) = KmTn (r) T„ (z) , (29) 

n=0 m,=0 
TV M 



p{r,z) = (^nniTn (t) (z) , (30) 

n=0 m=0 
N M 

e{r,z) = dnmTn (r) Tm {z) . (31) 

?i=0 m=0 

Expressions (p8|)-(|31|) are replaced into the equations ([T6|)-(p!9D and boundary condi- 
tions (^-(^. The N+l Gauss-Lobato points (r^- = cos(7r(l - j/iV)), j = 0, N) 
in the r axis and the M+1 Gauss-Lobato points {zj = cos(7r(l — j/M), j = 0, M) 
in the z axis are calculated. The previous equations are evaluated at these points ac- 
cording to the rules explained in Ref. in this way 4 (A^ + 1) (M + 1) equations 
are obtained with 4 (A^ + 1) (M + 1) unknowns. The system has not maximun rank 
because pressure is only determined up to a constant value. Since this value is not 
affecting to the other physical magnitudes, we replace the evaluation of the normal 
component of the momentum equations at (r^^^v = 1, -2j=4 = cos(7r(l — 4/M))) by a 
value for the pressure at this point, for instance p = 0. To solve the resulting nonlin- 
ear equatins, a Newton-like iterative method is used. In a first step the nonlinearity 
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is discounted and a solution is found by solving the linear system. It is corrected by 
small perturbation fields: Ur (r, z) , Uz (r, z) , p (r, z) and 9 (r, z) , 

<+i(r,z) = uiir,z)+Urir,z), (32) 

<+Mr,z) = uUr,z)+Uz{r,z), (33) 

(r, z) = (r, 2;) + p (r, 2;) , (34) 

e^+i(r,z) = {r, z) + Q {r, z) , (35) 

in such way that solutions at i + 1 step are obtained after solving Eqs. ([T6D -(p6|) 
linearized around the approach at step i. The considered criterion of convergence 
is that the difference between two consecutive approximations in norm should be 
smaller than 10~^. 

3.2 Results on the basic state 

Previous theoretical works [|T^, ID, dealing with this problem have approached 
the basic state solution as a parallel flow, where only the horizontal component 
of the velocity field exists. This approach requires the imposition of a constant 
temperature gradient over all the fluid layer and at the top boundary. In those 
descriptions this gradient becomes the main control parameter. In our model lateral 
effects are considered, so parallel flow is no longer valid and no restrictions are 
needed in the temperature boundary conditions which, as explained in the previous 
section, is the usual Newton law. This boundary condition is the origin of two 
control parameters related to temperature, i. e., AT^ and AT. This possibility. 



which had been already addressed in |21], is explored in detail in the next section. 



We have found two types of basic states. The first one, displayed in Fig. 2 
shows the isotherms obtained at i? = 1.25, 6* = 10 and AhT = AT = 6.4° C. The 
bottom profile is linear, as expected from the boundary condition, however at the 
top it approaches the constant ambient temperature, since the Biot number is large. 
This field is similar to that of the layer heated from below, i.e. roughly speaking 
hotter at the bottom. This feature is shared with the linear flow described in . 
Fig. 3 shows the velocity fields for the same conditions. They are formed by two co- 
rotative rolls perpendicular to the gradient. This result is similar to that obtained 
in experiments reported in Refs. |T^, 0, where co-rotative rolls perpendicular 



to the gradient are found as well. These results coincide with those obtained in Ref. 



21| for a rectangular geometry. 

The second kind of basic state is represented in Fig. 4 for B = 0.8, 6* = 2.5, 
AhT = 0.3° C, AT = 1.84° C. The most striking feature is that along the vertical 
axis there are changes in the sign of the temperature gradient. It seems that the 
Biot number is not dissipating heat effectively at the top and then heat is advected 
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by the velocity field. This temperature profile coincides with that of the return flow 
defined in Ref. |^|, where the role of the transport of heat by the velocity field 



as the origin of the instability is discussed. As we will demonstrate in the next 
section this flow becomes unstable through an oscillatory bifurcation, as it is also 
reported in |2D[. Crucial to the appearence of the two different flows described in 



this section, is not only the Biot number that dissipates heat at the top, but also 
the role of parameters AT and ATh which introduces heat into the system. 

Fig. 5 shows a transition between the two states we have just described. It 
is obtained at B = 0.3, 6* = 10, AhT = 10° C, AT = 20.41° C. At this stage the 
flow already becomes unstable through an oscillatory bifurcation. 

4 Linear stability of the basic state 

The stability of the basic state is studied by perturbing it with a vector field de- 
pending on the r, (f) and z coordinates, in a fully 3D analysis: 

Ur (r, 0, z) = ul (r, z) + Ur (r, z) e^-<^+^*, (36) 

(r, 0, z) = u\ (r, z) + (r, z) e^-*+^*, (37) 

(r, 0, z) = u\ (r, z) + (r, z) e^'^<^+"*, (38) 

e (r, 0, z) = (r, ^) + e (r, z) e''^*+^\ (39) 

p (r, 0, z) = / (r, z)+p (r, z) 6*°^*+^*. (40) 

Here the superscript b indicates the corresponding quantity in the basic state and 
the bar refers to the perturbation. We have considered Fourier modes expansions 
in the angular direction, because along it the boundary conditions are periodic. We 
replace the expressions (PBD-(P]D into the basic equations and after hnearizing the 



resulting system, we obtain the following eigenvalue problem (the bars have been 
dropped): 



AmUr — AdrP — G'^Ur — 2G'^imU(p 


= 0, 


(41) 


AmU(j> — Gimp + 2G'^imUr — G'^u^ 


= 0, 


(42) 


A^u^ - 2d^p + i?e 


= 0, 


(43) 


GUr + AdrUr + GimUff, + 2dzUz 


= 0, 


(44) 


A,^0 - UrAdrQ^ - u^Ade - 2u%e - 2uzdS^ 


= AG, 


(45) 



where Am = A'^d^ + GAdr — ra^G^ + 49f . The following boundary conditions for 
the pertubations are obtained, 

Uz = 2dzUr + MAdrQ = 2d,u^ + GMime = 2dzQ + BQ = 0, on ^ = 1,(46) 
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Ur = U(j, = Uz = Q = 0, on z = —1, (47) 
Ur = = M2 = 0, 9 = 0, on r = 1, (48) 
Ur = U(j, = Uz =, Q = 0, on r = —1, (49) 

together with 

AmMr — AdrP — 2G'^im.u^ = 0, on r = ±1, (50) 

AniM^ - 29^p + i?e = 0, on 2 = -1, (51) 

Gur + AdrUr + Gioiu^ + Id^Uz = 0, on 2; = 1. (52) 

4.1 Numerical method 

The eigenvalue problem is discretized with the Chebyshev collocation method used 
for the basic state. Now there is a new field, the angular velocity, which is expanded 
as follows, 

N M 

u^{r, z) = ^Y1 ^nmTn {r) (z) . (53) 

n=0 m=0 

In order to calculate the eigenfunctions and thresholds of the generalized eigenvalue 
problem the equations are posed at the collocation points according to the rules 
explained in [0, so that a total of 5(A^ + 1)(M + 1) algebraic equations are obtained 



with the same number of unknowns. If the coefficients of the unknowns which 
form the matrices A and B satisfy det(74 — XB) = 0, a nontrivial solution of the 
linear homogeneous system exists. This condition generates a dispersion relation 
A = \(m, R, M, B,u^,Q^,p^), equivalent to calculate directly the eigenvalues from 
the system AX = XBX, where X is the vector which contains the unknowns. If 
Re(A) < the basic state is stable while if Re(A) > the basic state becomes 
unstable. In the Appendix a study of the convergence of the thresholds calculated 
with this method is shown. 



4.2 Results on the stability analysis 

The basic solutions obtained in section 3, become unstable when the control pa- 
rameter, AT in this work, is increased. When AT is changed, AT^ which is not 
the control parameter in our choice, is fixed to a non zero value. In the explored 
parameter range two types of bifurcations take place: stationary and oscillatory. 
Figure 6 depicts the m dependence of Re(A) for B = 1.25, 6* = 10 and AT^ = AT 
at the threshold (AT^ = 6.4° C) and below it (AT < 6.4° C). The bifurcation is 
stationary since the imaginary part of the eigenvalue A is zero at the critical value 
mc = 20. Figure 7 is similar to Figure 6 but at B = 0.5, 6* = 10 and AT^ = 5° C 
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and ATc = 8.63° C. Now the bifurcation is oscillatory since the imaginary part of 
the eigenvalue is non zero at the critical value nic = 18. In order to understand how 
physical conditions affect the instability, Fig. 8 summarizes the influence of heat 
related parameters in the bifurcation. There, it is shown the critical value of ATc, 
as a function of the Biot number B, for several values of AT^. In this figure station- 
ary and oscillatory bifurcations are observed and it shows that having one or other 
transition depends on an equilibrium between these parameters. Travelling waves 
are favoured by systems storing a lot of heat, i.e., low values of the Biot number, and 
high values of A^T, which facilitates the transport of heat by the velocity field . 
We notice here the novelty of these results, where transitions between stationary 
rolls and oscillatory waves are due to changes in heat related parameters. This has 
not been reported in any experimental or theoretical work before, where transitions 



have been found due to variations in Prandtl number pO[ or in aspect ratio [1^, ^ 



In Fig. 8 it is also clear that thresholds decrease as B increases while the 
dependence of thresholds on ATh is not monotonous, it presents a minimum for 
AT/i = 2° C as figure 9 a) displays. In Fig. 9 b) it is shown how the critical 
wavenumber increases with AT/j, on the other hand it remains almost constant with 
B. 

The geometry and size of the box is also affecting the instability as it is 
noticed mainly in low aspect ratio systems. Although a detailed study of patterns 
in small containers is beyond the scope of this work, table I gives an insight into 
it. It shows the critical value of AT^ and the corresponding critical wavenumber m^ 
for some 6*. The wavenumber does not change monotonously with the aspect ratio, 
as also happens in small containers with uniform heating from below [1^, |11[] as 
well. 

The critical wavenumber for the stationary bifurcation shows that a 3D struc- 
ture appears which consists of radial rolls whose axes are parallel to the radial co- 
ordinate as Fig. 10 shows for B = 1.25 and 6* = 10 and ATc = 6.4°C. The growing 
perturbation along the transverse plane is plotted in Fig. 11, where a structure 



appears near the hot side as in the experiment reported in p3[. This is so because 



in this case, the instability has the same origin as a static fluid heated from below, 
so since the vertical temperature gradient is greater at this boundary, perturbation 
grows there first. On the other hand since the hot cylinder counteracts the vertical 
temperature gradient it inhibits the instability threshold. These results are in good 
agreement with those obtained in Ref. in similar conditions but in a cartesian 



geometry. In that article at a threshold ATc = 6.8°C, the critical wavenumber is 
kc = 2.4, which would correspond to 12 — 36 periods respectively in the inner and 
outer radii of an annulus as ours. This matches rather well with our critical thresh- 
old and wavenumber which is mc = 20, suggesting that, at least for this aspect 
ratio, geometry does not strongly affect the results. These results and discussions 
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of the stationary bifurcation are in good agreement with those reported in Refs. 



13, 18, 19 



We also obtain oscillatory bifurcations. For instance, there is one at 6* = 10 
[d = 2 mm), ATh = 5°C, B = 0.5 and ATc = 8.63°C, which corresponds to 
travelling waves or hydrothermal waves described in Refs. JI^ |T2|, |T^. Fig. 12 
depicts the eigenfunction along the transverse plane, where the tilted axis of the 
waves can be seen. This result is comparable to that reported in Ref. where 
for a cartesian container with = 20 mm, c? ~ 2 mm, hydrothermal waves of 
wavenumber fc^ ~ 1, emerge above a threshold. This would correspond to a pattern 
with 6 — 19 periods in the inner and outer radii of our annulus respectively, which 
is not far from our critical wavenumber m^ = 18. The angle of inclination of the 
waves in Fig. 10 is ~ Sir /A also close to that of 2.60 — 1.75 rad of Ref. [|l3 . 



5 Conclusions 



We have studied a BM lateral heating problem in a cylindrical annulus. The problem 
has been solved with a Chebyshev collocation method by keeping the original Navier- 
Stokes equations, where appropriate boundary conditions are required for pressure 
T0| . The scheme developed in Ref. |jTO| for stability problems has been extended 
to calculate a non trivial basic state. The procedure is confirmed to be reliable, 
effective and easy to implement. 

We have obtained two kinds of basic solutions. One type is formed by co- 

m 



rotative rolls such as those reported in experiments 12, ^ 



and in similar 



theoretical studies 21 



The second kind of solution has similar features to those of 
the return flow described in |pO|| . 

The linear stability analysis of these solutions shows stationary bifurcations 
to radial rolls and oscillatory bifurcations to hydrothermal waves, whose appearance 
has been proved to depend on their heat related parameters {B, AT^ and AT). This 
fact has not been adressed before, perhaps because experimentally these parame- 
ters are hard to manipulate. Properties of the growing patterns coincide with those 
reported in previous works. Stationary solutions appear in basic flows with temper- 
ature profiles close to the static fluid heated from below ^ . Their wavenumbers 
are in the range of experimental |l^, |13], [l^ and theoretical studies I^Tj. Also the 
structure in the z — r plane agrees with previous results ||2l|, ^ . Oscillatory solu- 
tions emerge from a basic state of the type of return flow as predicted in [^. Their 
structure is comparable to that obtained in experiments fl^ . 
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Appendix 



Convergence of the global numerical method 

To carry out a test on the convergence of the global numerical method (calculation of 
the basic state and subsequent linear stability analysis) we compare the differences 
in the thresholds of the differences of temperature (AT) to different orders of expan- 
sions for some values of the parameters involved in the problem: S* G [8.26, 10.52]. 
In table II the thresholds for these states varying the aspect ratio are shown for four 
consecutive expansions varying the number of polynomials taken in the r (A^) and 
z (M) coordinates. ^From the table it is seen that convergence is reached within a 
relative precision for ATc of 10~^. If M is increased the difference between succes- 
sive expansions is order 10~^ while increasing N it is order 10~^. Further increments 
of the order of the expansions give thresholds values oscillating around the quoted 
ones, so within that error, thresholds are convergent. The range of parameters has 
been selected in order to reach this convergence criterium. 
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Table captions 
Table I 

Critical AT and critical nic vs the aspect ratio 5*. The value of the horizontal 
temperature difference is AT^ = 0.4 and B — 1. 

Table II 

Critical temperature differences (°C) for different values of the aspect ratio at con- 
secutive orders in the expansion in Chebyshev polynomials {B — 1.25). 

Figure captions 

Figure 1 

Problem set up (a = 0.01 m, 5 = 0.02 m). 
Figure 2 

Isotherms of the basic state corresponding to values of the parameters R = 2228, M — 
51, S* = 10, B = 1.25 (AT = 6.4°C, AT,, = 6.4''C, d = 2 mm). 

Figure 3 

Velocity field at same conditions of Fig. 2. 
Figure 4 

Isotherms and velocity field of the basic state corresponding to values of the param- 
eters R = 40995, M = 59, 5* = 2.5, B = 0.8 (AT = 1.84°C, AT^ = O.S^C, d = 8 
mm). 

Figure 5 

Isotherms and velocity field of the basic state corresponding to values of the param- 
eters R = 7105, M = 163, 6* = 10, B = 0.3 (AT = 20.41''C, AT,, = 10°C, d = 2 
mm). 

Figure 6 

Maximum real part of the growth rate A as a function of m for basic state at M = 51, 
6* = 10, B = 1.25. The top line stands for the threshold condition R^ = 2228 (or 
ATc = 6.4°C), and its maximum determines the critical mg. The bottom line is for 
R < Re- Sohd lines correspond to branches with real eigenvalues while the dotted 
ones stand for complex eigenvalues. 

Figure 7 
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Maximum real part of the growth rate A as a function of m for a basic state at 
M = 69, 5* = 10, S = 0.5, at the threshold condition = 3004 (or AT^ = 8.63°C). 
The maximum determines the critical m^. The solid line correspond to a branch with 
real eigenvalues while the dotted one is for complex eigenvalues. 

Figure 8 

Critical AT vs the Biot number B for different values of AT^ [S* = 10). 
Figure 9 

Critical AT and m^ vs AT^ for 5* = 10 and B = 0.7. 
Figure 10 

Isotherms of the growing perturbation in the x — y plane at a stationary bifurcation 
point (AT = ATh = 6.4° C, 5* = 10, S = 1.25). 

Figure 11 

a) Isotherms of the growing perturbation in the r—z plane at a stationary bifurcation 
point (AT = ATh = 6.4° C, 6* = 10, B = 1.25). b) Velocity field of the growing 
perturbation in the r — z plane at same conditions. 

Figure 12 

a) Isotherms of the eigcnfunction in the x — y plane at an oscillatory bifurcation 
point {ATh = 5° C, AT = 8.63° C, S* = 10, B = 0.5). 

Figure 13 

a) Isotherms of the eigcnfunction in the r — z plane at the oscillatory bifurcation 
point considered in Fig. 12. b) Velocity field of the eigenfunction in the r — z plane 
at same conditions. 
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Table I 



s* 


AT, 




2.5 


3.91 


11 


10/3 


2.63 


1 


5 


3.28 


2 


10 


1.03 


10 



Table II 



s* 


13 X 9 


13 X 13 


25 X 13 


29 X 13 


10.52 


6.8711 


6.8519 


6.9551 


6.9229 


10.25 


6.5131 


6.5224 


6.6792 


6.6665 


10.00 


6.2223 


6.2298 


6.3717 


6.3938 


9.52 


5.8794 


5.9490 


5.9126 


5.9519 


9.09 


5.8415 


6.0807 


5.6626 


5.6868 


8.26 


5.7341 


5.9204 


5.6812 


5.6892 
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